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Abstract 



In this paper we give the distribution of the position of the particle in the asym- 
metric simple exclusion process (ASEP) with the alternating initial condition. 
That is, we find F(X m (t) < x) where X m (t) is the position of the particle at 
time t which was at m = 2k — 1, k € Z at t = 0. As in the ASEP with the step 
initial condition, there arises a new combinatorial identity for the alternating 
initial condition, and this identity relates the integrand to a determinantal form 
together with an extra product. 



1 Introduction 



The exclusion process is an interacting stochastic particle system on a countable set 
S. A particle at x € S chooses y € S with probability p(x, y) after a holding time 
exponentially distributed with parameter 1. If y is empty, the particle at x jumps 
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to y but if y is already occupied, then the particle remains at x, and the Poisson 
clock resumes. The detailed references on the construction of the model are Liggett 's 
books [21 E]- The asymmetric simple exclusion process (ASEP) is defined on S = Z 
by taking p(x, x + 1) = p and p(x, x — 1) = q for all i£Z, where p + q — 1. If p — 1, 
we call it the totally asymmetric simple exclusion process (TASEP). 

Schiitz [0] considered the system of N particles for the TASEP. There Schiitz ob- 
tained the probability that the system is in configuration • • • , xn} at time t given 
the initial configuration {y X) ■ ■ ■ , ?/Ar} at t = 0, and expressed the probability as an 
N x N determinant, and moreover, for p, q ^ 0, obtained the conditional probability 
for N = 1,2. Johansson |4] studied the TASEP with a special initial condition that 
all non-positive integer sites are occupied at t = 0, which is called the step initial 
condition. Johansson derived the probability that the iVth particle from the right- 
most moves at least M steps before time t, which describes the current. By using the 
technique used in [S] Rakos and Schiitz [S] obtained the same result as Johansson's [I] 
for the TASEP with the step initial condition. Rakos and Schiitz's method originates 
from the Bethe Ansatz while Johansson used a combinatorial argument. 

Tracy and Widom extended the previous results on the TASEP to the ASEP. The- 
orem 5.1 and Theorem 5.2 in [13] provide with the probability that the rath particle 
from the leftmost is at x G Z at time t when the system initially has finitely many 
particles. Corollary to Theorem 5.20 is a generalization of the Johansson's result and 
the Rakos and Schiitz's result to the ASEP and the subsequent Remark shows that 
how their result leads to the Johansson's result or the Rakos and Schiitz's result. The 
work for the deterministic step initial condition in (13] was generalized to the random 
step initial condition, so called, the step Bernoulli initial condition [15J. 

In the ASEP with the step initial condition, the random variable x m (t), the rath 
particle's position from the leftmost (or the rightmost depending on which side is 
occupied initially) at time t can be related to the current. If we assume the initial 
condition that all positive sites are occupied and other sites are empty at t — and 
P < Q, 

F(T(x,t) > ra) = F{x m (t) < x) 

where T(x, t) is the number of particles whose positions are less than equal to x at 
time t. The asymptotics on the fluctuation of x m (t) or T(x,t) are related to the 
Tracy- Widom distribution in random matrix theory and are well known for various 

1 Tracy and Widom assumed that initially all positive integer sites are occupied and other sites 
are empty. 
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situations HJ H El UM ■ 

Besides the step initial condition, the TASEP with the alternating initial condi- 
tion has been investigated [2JEJ. This initial condition assumes that all even sites are 
occupied and all odd sites are empty at t = or vice-versa. While the current fluc- 
tuation of the TASEP (more generally ASEP) with the step initial condition in the 
KPZ scaling limit is governed by the GUE Tracy- Widom distribution [Tj l%| [TT| [T5] . 
on the other hand, the current fluctuation of the TASEP with the alternating initial 
condition is related to the GOE Tracy- Widom distribution (2JEICE2]- But to the best 
of the author's knowledge nothing is known on the ASEP with the alternating initial 
condition although we expect the GOE statistics in this case. 

In this paper, we study the distribution of a particle's position in the ASEP with 
the alternating initial condition, which gives information on the current of the system. 
The asymptotic behavior of the current remains an unsolved problem. 

We denote by Y C Z the set of initial positions of particles and by Py the prob- 
ability of the ASEP with the initial condition Y. The main object we are interested 
in is a random variable X m (t), the position of a particle at time t whose initial 
position is m G Y and the main goal in this paper is to obtain Py(X m (t) < x) 
when Y = 2Z — 1 = {2i — 1 : i G Z}J^| To do so we will start with a finite set 

Y = {2i — 1 : i G Z, — N + 1 < i < N} and will consider the limiting case that 
iV — > oo. Let Y + = {1, 3, • • • , 2N - 1} and Y_ = {-1, -3, ■ ■ • , -2N + 1} so that 

Y = Y"_ |J Y + . Let us denote the position of the ith particle from the leftmost at time 
t by %i(t). Then, 

X m (t) = £[y_| + m±l(t), (1.1) 

so we can use some previous results on Xi(t) in the ASEP with finite Y. In Section 2 
we review the integral formula for Fyix^t) < x) when Y is finite and in Section 3 by 
using a new combinatorial identity that arises in the alternating initial condition we 
derive P 2 z_i(X m (t) < x) and give the one-sided version of it. The identity is given in 
Lemma 3.2 and our final results for P2z-i(A m (t) < x) and its one-sided version are 
given in (13.131) and (I3.15p . respectively. 

2 We consider 2Z — 1 instead of 2Z because of a technical reason. 
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2 Some known results 



The integral formulas for the probability of a particle's position for the ASEP with 
a finite initial condition Y were developed [T3l [T5| [16]. These integral formulas have 
different forms depending on the contours we choose. The integral formula of The- 
orem 5.2 in [13] which is over large contours can be used to derive the formula for 
the step initial condition that positive integers are occupied as shown in Corollary. 
Alternatively, if all negative integers are initially occupied and other sites are empty, 
the integral formula of Theorem 5.1 in [13] which is over small contours is needed. 
These are because we need geometric series which arise to be convergent. Hence, for 
the alternating initial condition that has infinitely many particles on both sides of 
any reference site, a recently developed integral formula over both large and small 
contours is required [16]. In this section we review the formula developed in [16] 
which is the starting point of this paper. 

First, let r = p/q and £ = £(k) = (fi, • ■ ■ ,6b + ,£-i, • • ■ ,£-k_) witl1 k = k + + k— 
We define 



efo) := 1+qti-l, /(6,fr):= A/* and I(x, := JJ /(& &) R ^_ , 

P + <KiSj s« i<c ^ 4 1 si 

Notice that I(x,£) depends on t but we omit it in the notation. Given two sets U 
and V of integers 

a(U, V) : = #{(«, v) : u G f7, v G V and u > u.} 
and recall the definition of the r-binomial coefficient, 



N 
n 



;i-r JV )(l-r 7V - 1 )---(l-r JV - n+1 ; 



(l-r)(l-r 2 )---(l-r«) 

We assume that FcZ and Y — Y_ (J 1 + is a finite set where Y"_ and Y" + are disjoint 
and all members of V + are greater than all members in F_. For S± C Y± we set 
|<S±| — k± and use positive indices for S+ and negative indices for S-. In other words, 
we set S- = {s-i, s_2, • ■ ■ , and 5+ = {si, sz, ■ ■ ■ , Sk + }. Then the distribution 

of x m (t), the mth particle's position from the leftmost particle at time t is given by 

PrOUt) < x) = E c m,s.,s + [ k ! /(' Oil- II < ( 2 - 2 ) 

|s±|=fe± 
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where 



-_ 1 Nm+|y_\s_| r m( ^- 1) -fc + m+ t 7(5 + ,y)+ f 7(y_,y-\5-)-m|y-|+ fc - (fc - +1) 



fc(fc-l) 

xq 2 



m 



k-l 

|y_ \ s_ 



and (C r ) is a circle with center zero and radius TZ (r). Here 1Z (r) is so large (small) 
that all the poles of Yli<j /(&>£?) i n I( x iO ne ms ide (outside) C-r. (C r ). The product 
in the integrand is over all positive and negative indices, and integrals over C-r are for 
variables of positive indices and integrals over C r are for variables of negative indices. 
If we set Y_ = 0, (I2.2p exactly becomes (3) in |16j . from which we are able to work on 
the (positive) one-sided step initial condition or the (positive) one-sided alternating 
initial condition. 

Additionally, in Section V of [TB], a variant of c m) s^,s+ was computed. It is the 
coefficient 



C m +\Y-\,S-,S-i 



<_l\m+k- T a(S + y + \S + )-a{Y-\S-,S-)+m(m-l)/2+ k + (k + +1) - 



(2.3) 



x *(fc-l)/2 



k-l 

m + k- — 1 



and will be used in the next section. 



3 Alternating initial condition on Z 

1 Distribution of a particle's position at time t 

In this section we derive P2z-i(A m (t) < x) as the limiting case iV — > oo of the formula 
fl2T2|) withF = {-2N+1, -2N+3, ■ ■ ■ ,2N-3,2N-1} and with X| y _| + m+i (t) instead 
of x m (t)i| Then S+ and S_, subsets of Y + and y_, respectively, may be written as 

S + = {2i 1 -l,2(i 1 +i 2 )-l,-.. ,2(i 1 + --- + i k+ )-l}, (i x ,... G N and A; + < N) 

and 

5_ = {-2JX + 1, -2(j 1+ j 2 )+l, ■ • • , -2(;i+- • -+J fc J+l}, (ji, • • • , jk- e N and < N). 
3 Recall (TO) . 
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Since we are working with X\ Y _\ + m+i(t) = X m (t), the coefficient for X m (t) becomes 

k- 1 



Xg k(k-l)/2 



m—1 



+ k_ 



by replacing m by ^ in (TOD . Noticing that Y~+ \ 5+) - a(Y_ \ S-, S-) in 

C\ Y _\+in±i,s-,s + anc ^ iLC* i n the integrand in (I2.2p depend on S±, we consider the 
sum 

s±cy±, i 

|S±|=fc± 

= ^ r <r(s + ,y + \s + ) Jj^.-«i . r - ff (y_\s_,s_) jJ^-*<_ 

S+CY+, i>0 S-CY-, i<0 

|5+|=fe + |S_|=fc_ 

Let us compute the first sum. First, 

£ 7 Si _ ^-2ii+l^-2ii-2i 2 +l . . . £~ 2il 2i fe+ +1 



i>0 



-2i 



fc4 



Now, observe that the number of points in Y + less than or equal to si = 2(ii + ■ ■ • + 
ij) — 1 is i\ + • • • + %i and the number of points in S + less than or equal to s/ is I. 
Hence, the number of pairs (si,y) with si > y where si G S + and y G Y + \ S + is 
«i + • — h i\ — I, and thus 

a(S + ,Y+\S+) = (i 1 -l) + (i 1 + i 2 -2) + --- + (i 1 + --- + i k+ -k+) 

= + (k+ — l)i2 H Vik+ ^ • 

Denoting the first sum by 



S+CY+, i>0 
|5+|=fc+ 



|S+|=fc+ 
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If we assume that Y + = {2i — 1 : % = 1,2, •••}, the sum implies geometric series, 
which converge because we choose large contours for variables with positive indices 
as shown in (12.21) . Hence, 



6 • • • ^fc+ 

((£i • • <k + ) 2 - t* + ) ((6 • • -6 + ) 2 - r^) ■ ■ ■ {e k+ - r) ' 
Likewise for negative indices, 

= ^2ii-1^2ii+2ia-l . . . ^ji+-+2j k _-l 

i<0 

= (£-1 ■ ■< /.• ) -a , • • ■e-*_) 3yi (e-a ■ • ■e-*-) 2 * ■ ■■*_£:• 

For cr(F_ \ £_) we set 

Si = —S-i, S + = —S- and Y + = — YL, 
and then it is easily seen that 

a(yi\5_,5_) = a(5 + ,r + \S + ) 

and so 

a(S + ,Y + \S + ) = C/i - 1) + Oi + - 2) + • • • + 0"i + • ' • + Jfc- - 

MA;_ + 1) 

Recalling that we choose small contours for negative indices so that geometric series 
for negative indices converge, one can obtain for y_ = {2i — 1 : i = 0, — 1, —2 • • • } 

|5_|=fc_ 



fc -( fc -+ 1 ) S-lS-2 ' ' - £-fc_ 

' ((r*- - (C-i • • -^J 2 ) (r^- 1 - (£-2 • • -£- fe J 2 ) ■ ■ • (r - £ fc _) 
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and let 



S-1S-2 ' ' ' 



((r*- - (e-i • • • e- fc _) 2 ) (r^- 1 - (£-2 • • • e_ fe j 2 ) •••(*-- e fc J 



Hence we obtained 



P 2Z _x(X m (t) < z) = J2 °m,k ± I k f k I{x,0<p4k-^-)v + {k + ,C+)Y[d^ (3.4) 



where 



Cm,k± 



, -,x™+l,l, ( m 2 -l , fc(fc + l) , , _m+l k ) Hk-1) 
( — 1) 2 + fe - r V 8 + 2 2 K +^g 2 



jfe-1 



m— 1 



(3.5) 



Remark 3.1. In Corollary in [T3] the identity (1.7) in [13J was used in obtaining 
Pz+(^m(^) = x) for the step initial condition. We give the P^+(x m (t) < x) to be 
compared with the case of the alternating initial condition in the later section^ 



T (k-m)(k-m+l)/2 

' ^ (1 + T) k( - k ~ l )k\ 

k>m 



k-l 
k — m 



(3.6) 



e, - ii 



n (1 4 )te _ T) ^-^ 



Moreover, the integrand in P^+(x m (t) < x) could be expressed as a determinantal 
form by using another identity (3) in [14|. It states that 

^z+(x m (t)<x) = ^c m>fciT / • j det(if £j))i<i ,j<hd£i ■ ■ • d£ k (3.7) 

— „ J Ct? J Ct? 



k>m 



where 



(3.8) 



p + q£? - £ 

and c m ,k, T is a constant depending on m, k, and r. Here we introduce the identity for 
the later use. 



Lemma 3.1. [Lfl 



det 



1 



n n-A)( 



&)(&-r) 



4 We slightly changed the form of (2) in [14] , 
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2 Symmetrization and Combinatorial identity 



As mentioned in Remark 3.1 a combinatorial identity was found to derive the integral 
formula of the distribution in case of the step initial condition and the integrand of 
the formula can be expressed as a determinant. This identity is associated with a 
special initial structure of the system, that is, the step initial condition. So we may 
expect to have a new identity associated with the alternating initial condition. In 
this subsection we find the new combinatorial identity^ and obtain an alternate form 
of (13.41) by using the identity. 

Lemma 3.2. Let r = - and p + q = 1. For k G N 

V + <lia(i)ia(j) - £<r(i) 1 



En 



T 



1 n 1+ ;:£ +&) n^ 



(1 + T)t (*-l)/2 11 T - fife 11 Q - T 



Proof. The equality clearly holds for k — 1. Denote the left hand side by Lfc(£i> " " > £fc) 
and the right hand side by Rk(£i, • • • >Cfe)) an d assume that the identity holds for 
k — 1, i.e, = Let cx(l) = I. We change the sum over all permutations 

in L k to the double sum over I = 1, 2, • - • , k and (cr(2), • ■ ■ , cr(/c)) G Sfc-i, that is, 
E CTe§fc = Ez fe =iE ((7 (2),.., CT (fe)) e s fc _ 1 - Observe that 



£ 2 ...p 2 n- k — t 2 t 2 ...t 2 T k 

W(1)W(2) Scr(fc) — r — S1^2 Sfc — T 

for all a G §fc and 

k . 

TT g + Qia(i)i<r{j) ~ _ 1 T P + ^(^(l) ~ TT P + g£q(i)£a(j) ~ ju(i) 

^(j) - ^(i) - r>7, ^oo _ ^« 

i,j=l,--- ,k i,j=2,---,k 



i,j=2,- ,k 



5 This identity was conjectured by Craig A. Tracy [TU]- 
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and 

t it t t n \" TT p + g^(0^(j) ~ 6ko 
•kfc-HSi, • ■ • , g-i,g+i, ■ ■ ■ ,qk) — 2-^ \ I 1 H7 

(a(2),..,a(fc))eS fe _ 1 i>i, ^ ^ 

i,j=2,- ,k 



X 



(C(2) ' ' ' C(fc) rfc )(C(3) ' ' ' C(fc) rk X ) " ' (£<r(fc) T ) 



Hence 



" " ' > 6fc) — 7272 12 a5^TT~ — <f * \ — ~Lk-i{€u ' " " > 6-1) 6+i) " " ' >60 



1 fc 

X^II — c II p. — "-^*-i(Ci) • • • ,6-1)6+1)* •• )C*;)) 



where the second equality comes from the induction hypothesis. Our goal is to show 
that 

i ^ 6 6 

2^ II — c % __ c — ~^-i(£i> ' " ' 6-i) 6+1, • • • ) Cfc) — -Rfc(^i) • • • > 60 



fc 



V~ > TT P + ^66 ~ 6 ^fc-ljCl) • • • ) 6-1; 6+1) " " " ) 6cJ _ C 2 C 2 t2 _ fc 

trjj 6-6 ^i,---,6) -^•••^ T - 

Recalling the form of Rk, what we want to show is 



A 



(l + r)*-i = - T) 1} ^-6 l + (3 ' 9) 

For some technical reasons we multiply by (p+<z6 2— 6)(l+ r— 26) both the numerator 
and the denominator of the right hand side, and then using p + q^f — 6 = (<?6 ~ p)(6 — 
1), (13. 9p becomes 

• • • il - r k _ n« g + g66 - 6 n 66 2fr - 1 - r 
(i + r)*- 1 "tT n^6-6 i i i i + r-6-6(6-i)(?6-p)' l ' j 
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Observe that 

P + q&z-tiTT T-zCi 2z-1-t 2q k - 1 U l tf 



(z) ■— TT - TT — 

for z — >■ oo. Thus, the integral of #(z) over a circle with sufficiently large radius 1Z is 
equal to the integral of — ' over the circle, so, 

[ g(z)dz = 2g fr " 1 JJe? = X) ***9{z)- 

Recall that r = p/q and p + q = 1. It is easy to see that Res z=1 g(z) = p k /q and 
Res z= £ g(z) = p k /q. Using p + q£f - f, = - p)(£, - 1) we obtain 

Res z=il g{z) = -r)_[J_ 



Finally, observing that 1 + r = we can see that Hes z=1+T _^g(z) = Res z= ^g(z), 
and thus (I3.9P is verified. This completes the proof. □ 

We use this new identity to symmetrize the integrand in ( 13. 4p . First, we write the 
integrand as 



i(x, o<p+(k + , z+)v-(k-, e-) = n f&> & n t^f^^ ^-(k-,t- 



n /(&,£i)nT37"^ + ' e 



X 



i<0,j>0 0<i<j i>0 

r -i-r -i-r Pe e(fi)t 

II f^^U-f^r^-^^ 

i<j<0 i<0 ^ l 



The first bracket is a function of variables with positive indices that is given by 
II /te,6)IlT^^^^) 

0<i<i i>0 ^* 

n 6 - 6 fr gv^; 

. ,. <4 o - 6 • • • e k+ - T k +ma ■■■a,- ■ ■ ■ (e k+ - 1 



X 
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and its symmetrization by the identity in Lemma 3.2 is 

c, tt nr g+Vfe)f 1 1 n 1+r ~^+^ am) 

Cfc +ll p + g ^ 6 1H i_ 6 q_ t ) 11 r _ a ^.iij 

where ^ 

Cjfc. 



£;+! (1 + r )fe + (fe+-i)/2- 
Now, (13. lip is written as 

r TT TT 1 n / g +1 e etf0t (6-r) N tt 1 + + j 

_ ~ Ck det r i i n r g^^^ i n l+r -& + & 

= c k+ det(K +{ ^)) 1<M<k+ n 

where 

gfc f = c fc+ -(-l)-^( W )-^- 1 J/ a g fci - 

and 

tlx e{i')t pi 



Likewise, for the second bracket, let & = £Zj • Then 

n /(&.<,) = n fel ' e - = n 

& - 6 



n 



0<i<J g+p6^-6 

and 

r-lr-3 ~ -(2fc_-l) 
/, * x Si ?2 



(r fe - -(£!•■■ 6_)- 2 ) - (6 ■ "GO-- 2 ) • • • 0" - (6J~ 2 ) 

fc_(fc_+i) fc_(fc_ + i) ~ -2 ~ -4 ~ -2fc_ ~ 

{r k - - (6 • ■ ■ £_)- a ) (t*- 1 - (6 ■ ■ ■ eL)~ a ) • • • (r - (£_)- 3 ) 

k_(fc_ + l) 

r — ^ — -(ei---e fc -) 

((6 • • ■ eL) 2 - Jr) (& ■ ■ ■ cL) 2 - ^=r) • • ■ ((gL) 2 - J) ' 
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Using Lemma 3.2 again, the symmetrization of 

k_(k_+l) 

n 



0<i<j 

is 



f Q + Viii ~ li ((6 " " " 6O 2 - ^r) ((6 ■ ■ ■ (k J 2 - t^t) • • ■ ((a- J 2 - J) 



Cfc_ 



11 I: r U II 



q. , ; 7 ' /'CO 6 j >0 — T 1 o<i<j T 1 ^0 

where 

fc_(fc-+i) 1 1 
Cfc - = r 2 (1 + r -i)fe-(fc--i)/2' 

and recalling that ^ = , we obtain the symmetrization of the second bracket 

(3.12) 

Using Lemma 3.1, (13.121) is expressed as 

cvdet ii 

i I) ^ 'J 

where 

c fe _ = c fc _-(-l)- fe -(p ? )- fe -( fc - 1 )/V- 

and 



Clx s(Z')t it'}- 1 

K-&0= -77, 7 -(Z'-T) U 



P + Qtt'-t v y (e')- 2 -T- r 

Now, we summarize our result. Recall the definition of 0) i n Section 2. Then 
we have 

Pzz-iPM*) < a;) (3.13) 

k ± >0 ^ Cr + i<0,j>0 i 



where 



k-<i,j<— 1 11 r -l _ t-ld~l 



= det^c^o)),™ n 1+r ~ ( £ +0 



i<i,i<fc+ li 
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and 



Cm,k± Ck-i- 



k,\kj. 



q 2 -l I fc+(fc+ + l) m+1 7. "\ fc(fc+l) 
~~8 1 2 ~~2~ K +) . 



q 2 x 



-k+(k+-l)/2 . -fc_(fe_-l)/2 



fc- 1 



m— 1 



3 The one-sided alternating initial condition 



Let us assume the initial condition Y' = {2n — ko '■ n G N} for a fixed k £ 1* . In 
this case we need to start with Theorem 5.2 @of O 

Pi-(^W<^) = X) c ™.*Z) rff(5,y) / ■"/ % fc '0llr^r"^ (3-14) 



ScY, 
\S\=k 



where 



fe(fe— 1)/2/ ijm^mfm-l)/ 



2— km 



k-1 
m — 1 



for convergence of geometric series. By using the same procedure as the subsection 
1 and the identity in Lemma 3.2 we can obtain Py/(x m (t) < x). We give the integral 
formula without the detailed derivation. 

_^ (fc-m)(fc~ro+l)/2 



F Y ,(x m (t) <x) = (-1) 



(i + r )fc(*-i)A;! 

k>m K ' 



k-1 

k — m 



x 



c n 



Zi - 6 



"+fco e te(? s ) 



n 



(3.15) 



This is to be compared with (13. 6 j) with the step initial condition and ( 13 . 13[) with the 
two-sided alternating initial condition. (" 13 . 15[) has a product term in the integrand 
compared with (13. 6p and is described by the integral only over large contours with 
the integrand G+(£) (when k = 1) compared with (I3.13p . 

Remark 3.2. One may consider more generalized periodic initial condition that sites 
in fcZ (k G N, k > 2) are occupied and all other sites are empty. The discrete TASEP 
with this initial condition was studied in [3]. We confirmed by using a computer 



3 This is an alternate form of the original one 
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that we do not have a combinatorial identity in a simple form when k > 3. That is, 
currently, we have two special deterministic initial conditions, that is, the step initial 
condition and the alternating initial condition, which are associated with combinato- 
rial identities. 
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Distribution of a particle's position in the ASEP 
with the alternating initial condition 

Eunghyun Lee 
September 28, 2010 

Abstract 



In this paper we give the distribution of the position of a particle in the asym- 
metric simple exclusion process (ASEP) with the alternating initial condition. 
That is, we find P(X m (t) < x) where X m (t) is the position of the particle at 
time t which was at m = 2k — 1, k € Z at t = 0. As in the ASEP with step 
initial condition, there arises a new combinatorial identity for the alternating 
initial condition, and this identity relates the integrand of the integral formula 
for ¥(X m (t) < x) to a determinantal form together with an extra product. 



1 Introduction 



The exclusion process is an interacting stochastic particle system on a countable set 
S. A particle at x e S chooses y 6 S with probability p(x, y) after a holding time 
exponentially distributed with parameter 1. If y is empty, the particle at x jumps 
to y but if y is already occupied, then the particle remains at x, and the Poisson 
clock resumes. The detailed references on the construction of the model are Liggett 's 
books [8j [9]. The asymmetric simple exclusion process (ASEP) is denned on S = Z 
by taking p(x, x + 1) = p and p(x, x — 1) = q for all x G Z, where p + q = 1. If p = 1, 
we call it the totally asymmetric simple exclusion process (TASEP). 
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Schiitz [12] considered the system of iV particles for the TASEP. There Schiitz 
obtained the probability that the system is in configuration {x\, ■ ■ ■ ,x^} at time t 
given the initial configuration {y\, ■ ■ ■ ,i/n} at t = 0, and expressed the probability 
as an N x N determinant, and moreover, for general ASEP, obtained the conditional 
probability for N = 1,2. Johansson [7\ studied the TASEP with a special initial 
condition that one half of the system is occupied and the other half of the system 
is empty at t — 0, which is called the step initial condition. Assuming the left half 
is occupied and the right half is empty at t = 0, Johansson derived the probability 
that the iVth particle from the rightmost moves at least M steps before time t, which 
describes the time-integrated current. By using the technique used in pjj Rakos 
and Schiitz [TT] obtained the same result as Johansson [7] for the TASEP with step 
initial condition. Rakos and Schiitz's method originates from the Bethe Ansatz while 
Johansson used a combinatorial argument. 

Tracy and Widom extended the previous results on the TASEP to the ASEP. 
Theorem 5.1 and Theorem 5.2 in [T6] provide the probability that the mth particle 
from the leftmost is at x G Z at time t when the system initially has finitely many 
particles. Corollary to Theorem 5.2 is a generalization of the Johansson's result and 
the Rakos and Schiitz's result to the ASEP and the subsequent Remark shows that 
how their result leads to the Johansson's result or the Rakos and Schiitz's result. 
The work for the deterministic step initial condition in [16] was recently generalized 
to step Bernoulli initial condition [18J that assumes that at t — each site in Z + is 
occupied with probability p with < p < 1 independently of the others and all other 
sites are empty. 

In the ASEP with step initial condition that the right half of the system is occupied 
and the left half is empty at t — 0, the random variable x m (t), the mth particle's 
position from the leftmost at time t can be related to the time-integrated current. 
Assuming p < q, 

F(T(x,t) > m) = P(x m (t) < x) 

where T(x, t) is the number of particles whose positions are less than equal to x at 
time t. That the asymptotics on the fluctuation of x m (t) or T(x, t) are related to the 
GUE Tracy- Widom distribution in random matrix theory and that fluctuations are 
in the regime of the t 1 ^ 3 scale are well known for various situations [21 0, dH HE] ■ 

Besides the step initial condition, the TASEP with the alternating initial con- 
dition also has been investigated. This initial condition assumes that all even sites 
are occupied and all odd sites are empty at t = or vice-versa. While the current 
fluctuation of the TASEP (more generally ASEP) with step initial condition in the 
long limit is governed by the GUE Tracy- Widom distribution [21 El [HI [TH [T5] , on the 
other hand, the current fluctuation of the TASEP with the alternating initial condi- 
tion is related to the GOE Tracy-Widom distribution [3-6,10,15]. The appearance of 
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the GOE Tracy- Widom distribution in growth models goes back to [1] by Baik and 
Rains and a breakthrough in the TASEP was led by Sasamoto [13]. But to the best 
of the author's knowledge nothing is known on the ASEP with the alternating initial 
condition although we expect the GOE statistics in this case. 

In this paper, we study the distribution of a particle's position in the ASEP with 
the alternating initial condition, which gives information on the current of the system. 
The asymptotic behavior of the current remains as a problem for the future. 

We denote by Y C Z the set of initial positions of particles and by Py the prob- 
ability of the ASEP with the initial condition Y. The main object we are interested 
in is a random variable X m (t), the position of a particle at time t whose initial 
position is m G Y and the main goal in this paper is to obtain Py(X m (t) < x) 
when Y = 2Z — 1 = {2i — 1 : % G Z}. To do so we will start with a finite set 

Y = {2i — 1 : i G Z, —N + 1 < i < N} and will consider the limiting case that 
N -> oo. Let Y + = {1, 3, • • • , 2N - 1} and F_ = {-1, -3, • ■ • , -2N + 1} so that 

Y = Y- [J Y + . Let us denote the position of the zth particle from the leftmost at time 
t by Xiit). Then, 

X m (t) = X\y_\ + m±l(t), (1.1) 

so we can use some previous results on Xi(t) in the ASEP with finite Y. In Section 2 
we review the integral formula for Py(xj(t) < x) when Y is finite and in Section 3 by 
using a new combinatorial identity that arises in the alternating initial condition we 
derive P 2 z-i(A m (t) < x) and give the one-sided version of it. The identity is given in 
Lemma 3.2 and our final results for P2z-i(A m (i) < x) and its one-sided version are 
given in (13.131) and (I3.15p . respectively. 



2 Some known results 



The integral formulas for the probability of a particle's position for the ASEP with 
a finite initial condition Y were developed [T6l [T8| [T9] . These integral formulas have 
different forms depending on the contours we choose. The integral formula of Theorem 
5.2 in [16] which is over large contours can be used to derive the formula for step initial 
condition that positive integers are occupied as shown in Corollary. Alternatively, if all 
negative integers are initially occupied and other sites are empty, the integral formula 
of Theorem 5.1 in [16] which is over small contours is needed. These are because 
we need geometric series which arise to be convergent. Hence, for the alternating 
initial condition that has infinitely many particles on both sides of any reference site, 
a recently developed integral formula over both large and small contours is required 
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j. In this section we review the formula developed in [19J which is the starting 
point of this paper. 

First, let r = p/q and £ = £(k) = (£1, • ■ ■ ,£k+,£-i, ■ • • ,£-k_) with k = k + + k— 
We define 

and 
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Notice that depends on t but we omit it in the notation. Given two sets U 

and V of integers 

a(U, V) : = #{(«, v) : w G U, v G V and u > u.} 
and recall the definition of the r-binomial coefficient, 

(l-r JV )(l-r iV - 1 )---(l-r Ar -' t + 1 ) 



N 
n 



(l-r)(l-r 2 )---(l-r") 



We assume that FcZ and Y = Y_ (J F + is a finite set where Y"_ and V + are disjoint 
and all members of y + are greater than all members in F_. For S± C Y± we set 
|<S±| = k± and use positive indices for S + and negative indices for 5*_. In other words, 
we set S- = {s-i, s_2, • • • , s~k_} and 5+ = {si, S2, ■ ■ ■ , $k + }. Then the distribution 
of x m (t), the mth particle's position from the leftmost particle at time t is given by 



where 



Cm,S-,S+ 



k±>0 S±CY±, 
\S±\=k± 



, 1 \m+|y-\5-| T ro( ^~ 1) -fc+m+ g (5 + ,y)+ CT (r-,y-\5-)-m|y-|+ fc - (fc - +1) 



fc(fc-l) 
xg 2 



jfe-1 

- 1 



and Cr. (C r ) is a circle with center zero and radius 7£ (r). Here 7£ (r) is so large (small) 
that all the poles of Yli<j * n lie inside (outside) Cr (C r ). The product 

in the integrand is over all positive and negative indices, and integrals over Cr are for 
variables of positive indices and integrals over C r are for variables of negative indices. 
If we set Yl = 0, (12. 2p exactly becomes (3) in [19] . from which we are able to work on 
the (positive) one-sided step initial condition or the (positive) one-sided alternating 
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initial condition. 

Additionally, in Section V of [TJJ], a variant of c mt s_,s+ was computed. It is the 
coefficient 



Cm+\Y-\,S-,S+ 



'_ 1 } m +k- r *(S + XAS+)-v(Y-\S-,S-)+ m( ™- 1) + k + ik + +1) -mk + 



x fc(*-l)/2 



jfe-1 

m + A;_ — 1 



(2.3) 



and will be used in the next section. 



3 Alternating initial condition on Z 



1 Distribution of a particle's position at time t 

In this section we derive $ > 2Z-i{X m (t) < x) as the limiting case N — > oo of the formula 
d22D with Y = {-2N+1, -2N + 3, ■ ■ • , 2iV-3, 2iV-l} and with a7 |y _ | , m +i (t) instead 
of x m {t). Then S + and S_, subsets of Y + and Y_, respectively, may be written as 

S + = {2i x - 1, 2{h + i 2 ) - 1, • • • , 2(«! + • • • + i k+ ) - 1}, (zi, • • • , i k+ E N, k + < N) 

and 

5_ = {-2j 1 + 1,-2(j 1 +j 2 ) + 1,--- ,-2(j 1 + ---+j fe _) + l}, {j 1 ,---J k _ EN, k^<N). 
Since we are working with X\ Y _\ + m±i(t) = X m (t), the coefficient for X m (t) becomes 

r - ( l)^+k-Jv(S + XAS + )-v(Y-\S-,S-)+^+ k+{k + +1) -^k + ) 

C\y_\ + hi±1,S-,S + - 2 rV " 

k- 1 



Xg fc(fe-l)/2 



m— 1 



2 

by replacing m by ^ in (Q • Noticing that cx(S + , F+ \ 5+) - \ SL, 5_) in 
C|y_| | m+i ) 5_ ) 5 + and riiC 7 ^ i n ^ ne integrand in (12. 2p depend on S±, we consider the 
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sum 



7 . -(s +I y + \s + )- ff (y_\5- I 5-) JJ^- 



\S±\=k± 



|S+|=fc+ |S_|=fc_ 
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Let us compute the first sum. First 

L-2i 2 + l 



i>0 



Now, observe that the number of points in Y + less than or equal to si = 2{i\ + • • ■ + 
ii) — 1 is i\ + • ■ ■ + i\ and the number of points in S + less than or equal to si is I. 
Hence, the number of pairs (si,y) with Si > y where Si G S+ and y G Y + \ S + is 
h H h %i — I, and thus 

a(S + ,Y+\S+) = (i 1 -l) + (i 1 + i 2 -2) + --- + (i l + --- + i k+ -k+) 

= + («+ — 1)«2 H • 

Denoting the first sum by f + (k + ,!; + ), 

Mk + ,t + ) = £ ^ (s+,ms+) I]C :S ' 

5+CY+, i>0 
|5+|=fc+ 

|S+|=fc+ 

If we assume that Y + = {2i — 1 : i = 1,2, • ■ ■ }, the sum implies geometric series, 
which converge because we choose large contours for variables with positive indices 
as shown in (12.21) . Hence, 

£l • • • £fc+ 

((£i • • " ((6 • • "6 + ) 2 - r^- 1 ) • • • {e k+ - r) ■ 

Likewise for negative indices, 

£ 7 Si _ ^2ii-1^2Ji+2i2-l . . . r' 2 y 2 -?'* 1 

= • • - e-fcj-^-i • • -e- fe j 2ji (e-2 ■ ■ - e-fcj^ ■ • -dt. 

For cr(Y_ \ £_) we set 

5; = — s-i, S + = —S- and Y + = — YL, 
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and then it is easily seen that 

a(Y_\S_,S_) = a(S + ,Y + \S + ) 

and so 

<r(S + ,Y + \S + ) = (j 1 -l) + (j 1 +j 2 -2) + --- + (j 1 + ---+j k _-k-) 

i ■ n ix- ■ + 1) 
= fc_j! + (A;_-l)j 2 + ---+j fe _ ^ '-. 

Recalling that we choose small contours for negative indices so that geometric series 
for negative indices converge, one can obtain for K_ = {2i — 1 : i = 0, — 1, —2 • • • } 



E 



T 



- CT (y_\5_,5_) 



n<r* 



5_cy_, 

|5-|=fc- 

fc_ <fc_+i) 
r 2 



fe_(fe_+i) 

r 2 



j<0 



E 



jii- jfc_=i 



(^_i • • -e-kj 2 vv(£-2 • • •£-*-) V v* 



fc_-i 



S-1S-2 ' ' ' S-fc_ 



- (£-x ■ ■ • ^J 2 ) (r^- 1 - (£_ 2 • • ■ £- fe _) 2 ) • • • (r - £j 



and let 



S-1S-2 ' ' ' 



((r fe - - (£-x • ■ ■ Z-k-) 2 ) (r*- 1 - (e- 2 ■ • • t-*_) 2 ) ' ' ' (r - £ J ' 
Hence we obtained 

W2z-i(X m (t)<x)=J2 c m,k±[ h [„ I(x,0^~(k-,UMK^ + )Hd^ (3.4) 

fc ± >0 JCn h+ Jc r j 



where 



m + l 



c m , fc± = (-l)^" +fe -r 



fc-1 



TO— 1 



(3.5) 



Remark 3.1. In Corollary in [16] the identity (1.7) in [16] was used in obtaining 
^z+ixmit) = x) for step initial condition. We give the F%+(x m (t) < x) to be compared 
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with the case of the alternating initial condition in the later section. 



(fc-m)(fc-m+l)/2 

F z+ (x m (t)<x) = (~l) m E ( 1+T w*-i )jb , 



k-l 

k — m 



(3.6) 



x 



& - 6 



ria 



Moreover, the integrand in F^+(x m (t) < x) could be expressed as a determinantal 
form by using another identity (3) in [T7]. It states that 

{x m {t)<x) = '^2 C m,k,T / • ■ / det(i^(^j, £j))l<ij<k 

( 3 - r ) 

fc>m i 



where 



x e{£)t 



(3.f 



and c mj fc jT is a constant depending on m, k, and r. Here we introduce the identity for 
the later use. 

Lemma 3.1. J7?|/ 
det 



1 



(-1)*(P9)^(?-*II 



-& 11 1-6. ( 



6) (6 



r 



2 Symmetrization and Combinatorial identity 



As mentioned in Remark 3.1a, combinatorial identity was found to derive the integral 
formula of the distribution in case of step initial condition and the integrand of the 
formula can be expressed as a determinant. This identity is associated with a special 
initial structure of the system, that is, the step initial condition. So we may expect to 
have a new identity associated with the alternating initial condition. In this subsection 
we find the new combinatorial identity^] and obtain an alternate form of f)3.4p by using 
the identity. 

1 This identity was conjectured by Craig A. Tracy through private communication. 



s 



Lemma 3.2. Let r = 2 and p + q = I. For k G N 

P + <?£a(i)6r(j) - 



En 



x 



(C(1)C(2) ' ' ' ^l(k) rfc )(^a(2)C(3) ' ' ' %l(k) X ) " ' 

1 TT 1 + T - & + 



II"- f 



(l +r) ^-l)/2ll r _^. 11 £2 



Proof. The equality clearly holds for /c = 1. Denote the left hand side by Lfc(£i, • • • , 
and the right hand side by Rk{£i, • • • ,£k), an d assume that the identity holds for 
k — 1, i.e, = Rk-i- Let <r(l) = Z. We change the sum over all permutations 

in Lfc to the double sum over / = 1, 2, • • • , and (cr(2), • • • , cr(/c)) G S^-i, that is, 
E, eSfc = Ef=iE( ff (2),..., ff (fc)) e s fc _ 1 - Observe that 

Ser(l)Scr(2) Ser(fc) r ~ S1S2 Sfc r 

for all a G Sfe and 

TT P + qC(»)6r(j) - 6r(i) 
TT P + - 6t(») TT P + Q^(i)^(j) ~ 

f = ' ^(i) - &(*) 4i ^ ~ £ ct « 

TT P + ~ 6 TT P + g&(i)&r(j) ~ &r(i) 



and 



1 

X 



(C(2) ' ' ' rfc )(C(3) ' ' ' C(fc) ^ " ' (^a(fe) T ) 
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Hence 

^fc(Ci) ■ • • > €k) — T2P2 — 72~z — fcX^II — c _ c — • • ■ ) 6-1; 6+1; " • >£*) 

Slfe""?* r l=1 ^, ?i Si 



A; 



~~ 7272 72 _ ^2 II ~ — — ' ' ' ) 6-1) 6+1) ■ • • >6s)) 

?i?2"--?fc T l=1 i+l Si Si 

where the second equality comes from the induction hypothesis. Our goal is to show 
that 

7172 — 7tz_ — k II — c _ <r — ^fe-i(6) ■ ■ ■ > 6-i) 6+i) • ■ ■ ) 60 = -Rfc(6) • ■ ■ ) 6c) 
S1S2 ■■■Sfe t- l=l & 

i.e., 

fc 

, 1 , 6 _ 6 -Rfc(6> • • • >&) 

Recalling the form of what we want to show is 

6 2 4l ■■■(t-T k A, 2 , rr p + g66 - 6 r - 66 

(l + r)*-i = - T) 1} ^-6 l + (3 ' 9) 

For some technical reasons we multiply by (p+<z6 2— 6)(l+ r— 26) both the numerator 
and the denominator of the right hand side, and then using p + q^f — 6 — (<?6 — p)(6 ~~ 
1), (13. 9p becomes 

sl£ 2 2 • • ■ 4 2 - r fc _ A n g + g66 6 tt r - 66 26-1-t 
(i + r)*- 1 "tT n^6-6 A l i i + r-a-6(6-i)(g6-p)' l ' ] 

Observe that 

p + ?6*-6tt r -<*6 2*-i-r 2g fc - 1 n I tf 



\^ 77 v + g66 ~ 6 -Rfc-i(6) • • • ) 6-i) 6+i) • • • ) 6c) _ C 2 C 2 & _ k 

2^11 c, _ t. P,YA ... t,^ -SiS2---S fe r . 



=11 ,_ 6 LLr 



2-6 AA 1 + r - z - 6 (z - l)(qz - p) 



for z 00. Thus, the integral of g(z) over a circle with sufficiently large radius 1Z is 

2„fe-l FT £2 

equal to the integral of 1 - !J — L over the circle, so, 
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Recall that r = p/q and p + q = 1. It is easy to see that Res z=1 g(z) = p k /q and 
Res 2=£ g(z) = p k /q. Using p + qtf - 6 = - p)(6 - 1) we obtain 



Res z= t l9 (z) = (Cf -r)[| 



p + g&fr - 6 t - 66 

6-6 1 + T - 6 - 6 



Finally, observing that 1 + r = we can see that Res z= i +T ^ l g(z) = Hes z= ^g(z), 
and thus ( 13. 9ft is verified. This completes the proof. □ □ 



We use this new identity to symmetrize the integrand in ( 13 ,4p . First, we write the 
integrand as 



i<j ^ 



i<0,j>0 0<i<jr' i>0 ^* 

r -i-r -i-r Pe e(?i)i 



I J " ' " 



The first bracket is a function of variables with positive indices that is given by 

0<i<j i>0 1 ^ (K^j P + * i>0 1 ^ 



n 



v + - 6 i 



6 - 6 m • • • e k+ - r^mn ■ ■ ■ e k+ - ■ ■ ■ (e k+ - r) 

and its symmetrization by the identity in Lemma 3.2 is 

0. '^'^ ^ i>0 1 * T 0<i<j ' ^< 

where 

C'k 



k + \ (i + r )Mfef-i)/2- 
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Now, (13. lip is written as 

TT 1 TT /£f +1 e e(6)t te-r)\ ^ 1 + r - + 



Cfc, 



TT TT * TT ^ VSi ~ ' M TT 

11 r,4_ rt A.t. _A.lln -AUt.-T^ 11 V P 2 - -r / H 



- Cfc+de H P+g ^-^A^ + iU c 2 -- >ll r-te, 

=~c k+ *et(K + &,Q) 1<hj<k+ n 1+ ;: ( 5* +ej) 

0<i<j ^ 3 

where 

c, + = c k+ ■ (-l)-*-(pq)-*+^-Wq* 

and 

tlx p e(?)t H 



Likewise, for the second bracket, let £j = • Then 

n /(&&) = n ^-'j-i n 6-6 



i<i<o o<j<i P + 9& £j _ 6 o<i<« P&j + 9 & 

= TT ~ & 



and 



r-lr-3 ~ -(2*_-l) 

(t*- - (6 • ■ ■ eL)- 2 ) (r^- 1 - (6 ■ ■ -£_)- 2 ) ■ ■ • (r - fe> 2 ) 

fc_(fc_+l) fc_(fc_ + i) ~ _2 ~ -4 ~ -2fc_ ~ 

^ 2 2 -ii 6 •••6- •(^•••^-) 

(r fc - - (6 ■ ■ • 60~ 2 ) (r*- 1 - (6 " " " CL)- 2 ) ■ ■ ■ (t - fej- 2 ) 

((6 " • " 6 J 2 " ^) ((6 • ■ " &_) 2 - ^=r) ■ • ■ (feJ 2 - i) ' 
Using Lemma 3.2 again, the symmetrization of 

k_(k_+l) 

TT k ~ & r ~ 2 -(6 •••6-) 
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IS 



n V/' c- n^-^— ) n i+t ~'- (£ ' +6 



„. , .. , 7 ' /'Cs.i 6 j >0 & — r 1 o<i<j r 1 
where 

fc-(fc-H-l) 1 1 



£;_! (1 + r - !)*-(*-- 1)/2' 
and recalling that £j = £ ^ , we obtain the symmetrization of the second bracket 

r n ^ rrf S^a. 1 a n 1 + 1/r ~ (cl + ^ 1) 



(3.12) 

Using Lemma 3.1, (13.121) is expressed as 

cvdet^ ii 

j<i<0 ^* ^3 

where 

5 fc _=c fc _-(-l)- fc -( M )- fc -( fe - 1 )/V- 

and 

fe E ^« (f')-i 



Now, we summarize our result. Recall the definition of £j) in Section 2. Then 
we have 

P 2Z _x(X m (t) < x) (3.13) 



fc±>0 "' C ^'' + i<0,j>0 i 



where 



and 



fe-<ij<— l 11 r -i _ t-it- 1 

i<j<o 

G+K+ ) = det(if + (6,&)) 1<i , <l+ n 1+ ;: ( 5; +fe) 

0<i<j ^3 



_ _ , „ v '"+1 1 f ^ 2 ~! I M fc + +1 > w+lt., ^ k(k+l) 

Cm,k ± C k+ C k . = (-1)~ + • fc , fc j 8 + 2 2 2 X 



-fc+(fc+-l)/2 . -*_(*_ -l)/2 



k-l 
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3 The one-sided alternating initial condition 



Let us assume the initial condition Y' = {2n — ko : n G N} for a fixed k$ G Z. In this 
case we start with (3) in [19], 



Py(xm(t) < a;) 



fc>l sty, i i 



sty 

\S\=k 



where 



C m ,k = q Hk-D/2^m T m(m-l)/2- 



km 



k-1 
m — 1 



for convergence of geometric series. By using the same procedure as the subsection 
1 and the identity in Lemma 3.2 we can obtain Py/(:c m (t) < x). We give the integral 
formula without the detailed derivation. 



F Y ,(x m (t)<x) = (-irJ2-R 



(k-m)(k-m+l) /2 



X 



n 



k>m 

U - k 



(1 + r )fc(*-i)jfe! 



k-1 

k — in 



^+fco p te(?i) 



i<j 



(3.15) 



This is to be compared with (13 .6p with step initial condition and ( 13 . 13[) with the 
two-sided alternating initial condition. (13. 15ft has a product term in the integrand 
compared with f 1 3 . 6 1) and is described by the integral only over large contours with 
the integrand (when k = 1) compared with (I3.13p . 

Remark 3.2. One may consider more generalized periodic initial condition that sites 
in fcZ (k G N, k > 2) are occupied and all other sites are empty. The discrete TASEP 
with this initial condition was studied in jl]. We confirmed by using a computer 
that we do not have a combinatorial identity in a simple form when k > 3. That is, 
currently, we have two special deterministic initial conditions, that is, the step initial 
condition and the alternating initial condition, which are associated with combinato- 
rial identities. 
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